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In this letter we apply dynamical system methods to study all evolutional paths admissible 
for all initial conditions of the FRW cosmological model with a non-minimally coupled to 
gravity scalar field and a barotropic fluid. We choose "energy variables" as phase variables. 
We reduce dynamics to a 3-dimensional dynamical system for an arbitrary potential of the 
scalar field in the phase space variables {K(j)/^/6II, kVV /VSH, K(f>/ \/6)- After postulating the 
potential parameter F as a function of A (defined as —V'/V) wc reduce whole dynamics to a 
3-dimensional dynamical system and study evolutional paths leading to current accelerating 
expansion. If we restrict the form of the potential then we will obtain a 2-dimensional 
dynamical system. We use the dynamical system approach to find a new generic quintessence 
scenario of approaching to the de Sitter attractor which appears only for the case of non- 
vanishing coupling constant. 

PACS numbers: 98.80.-k, 95.36.-hx 



I. INTRODUCTION 

At present a scalar field (j) plays a very important role in cosmological investigations. The 
discovery of cosmic acceleration 11,0] is a motivation to study dynamical models of dark energy 
which can be treated as some alternatives to the cosmological constant (see for review). In this 
context the simplest dynamical models involving the scalar field (j) (the quintessence idea 0, Q]) 
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minimally coupled to gravity with the potential V{(j)) are considered to model a time dependent 
equation of the state parameter w = p^/p^, where = (j)'^/2 — V{(j)), = 4)^/2 + V{(l)) and dot 
denotes differentiation with respect to the cosmic time t. 

On the other hand the unknown nature of dark energy expressed in terms of w can be constrained 
by a variety of astronomical observations. Recently the WMAP experiment has published its five- 



y, Q] and new supernovae datasets of Union compilation 
has also been published [8|]. It is interesting that some inflationary scenarios, i.e. forms of potential 
function of the scalar field can be rejected 



year data and polarization power spectra 



y, Q]. Furthermore, the single slow-rolling scalar field 
with potential V{(j)) oc 0^ is well within a 2a confidence level region, whereas another scalar 
field potential V{(j)) oc 0^ has been excluded on more than a 2a confidence level 3, y, [lo| . 
While the simplest candidate for dark energy seems to be a positive cosmological constant, i.e. 



n[Q, Q, 



such an 



the LambdaCDM model is favored by a Bayesian model selection method 
explanation of the accelerated expansion of the universe as associated with the vacuum energy 
meets the so-called fine tuning problem [l^ and it also suffers from the coincidence problem 
[isl . Therefore various alternative routes have been proposed like phantom or extended 



quintessence 



201, 
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22l | . While the minimally coupled scalar field endowed with a quadratic 
potential function has a strong motivation in inflationary models its generalizations with a simple 
non-minimal coupling term ^RcjP' have been studied 23|] in the context of the origin of the canonical 



inflaton field itself. The physical motivation to study of non-minimally coupled scalar field could 



be possible application of this mode 
has a long history (see for example 



s to inflationary cosmolo gy o r to the present dark energy and 
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281, 



291, 
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3d, 



33|, 



341, 
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36|, 
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38|, 



In the present paper, our aim is to perform a complete study of the global dynamics, attractor 
properties and stability of trajectories of both canonical and phantom scalar fields in the framework 
of dynamical systems methods. Recently the dynamical system methods was used to find attractor 
properties of phantom scalar field 40(]. Such an analysis was restricted to the case of the simplest 



ill). Our 



exponential or quadratic potential functions and minimally coupled scalar field (see also 
study will include a non-minimal coupling of the scalar field to the gravity in the form of the term 
and we will not consider any a priori form of the potential of the scalar field. We take as phase 



variables the same variables used in 



42 



43( 1 in the context of minimally coupled canonical scalar 



field. The first attempt to study the minimally coupled scalar field with a quadratic potential 
cosmology in these variables was made by Belinsky et al. 4J]. In this parameterization of the 



phase space the normalized variables of the scalar field and its first cosmic time derivative are 
used. The third new variable is related to the Hubble function and the resulting system is a 
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3— dimensional autonomous one. Earlier, a minimally coupled massive scalar field in the closed 



451] where analytical expression for the matter dominated 



FRW model was studied by Starobinsky 
and quasi-deSitter stages was derived. 

In our approach extended to the case of a non-minimally coupled scalar field we parameterize the 
phase space using the energy variables determined by the Friedmann energy equation. The choice 
of two variables is the same like in the Belinsky approach, but we complete the set of variables by 
a normalized variable related to the so-called roll parameter A oc V' IV , ' = d/d(h. Therefore our 
dynamical analysis is the generalization of the studies [4g, |47|] to the case of non- vanishing coupling 
constant ^. 

II. DYNAMICS OF NON-MINIMALLY COUPLED SCALAR FIELD COSMOLOGY IN 

THE ENERGY PHASE SPACE 

In this section the dynamical system methods are used to study the dynamics of the scalar 
field non-minimally coupled to gravity in the Robertson- Walker geometry. To make the discussion 
less complex, we ignore any coupling of the scalar field to matter and assume that the universe is 
spatially flat. We choose as phase space variables the so-called "energy variables" introduced in 



42l | which in the case of the minimal coupling (.^ = 0) assume the form 



where V{(j)) is a potential function of the scalar field (j), H = d/a is the Hubble parameter and a is 
the scale factor. These variables determine the ratio of kinetic and potential energy of the scalar 
field to the total energy 

where = ^H"^ for the flat model filled with the scalar field only. 

If we additionally postulate the presence of matter which satisfies the barotropic equation of 
state, then + pm = ^H'^j where pm = WmPm and Wm is constant and energy conservation 
requires 

^ + ^ = x^ + 2/^ + a. = l, (3) 

where VLm is the matter density parameter. Because of ([3]) the trajectories in the phase space will 
always be located within the unit circle in the phase space. If we restrict ourselves to positive and 
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monotonically decreasing potential functions V{(f)) we can consider only the trajectories located in 
one quadrant of the energy phase space. 

Above considerations can be simply generalized to the case of the non-minimal coupling of the 
scalar field (both canonical and phantom) to the curvature. In this case the action for the scalar 
field and gravity assumes the following form 

^=\J d'^^^ (^^^-^ {9'"'d^(t^d.(f> + ^R^^) - 2U{(I>)^ (4) 
where = SttG, £ = +1, — 1 corresponds to canonical and phantom scalar fields respectively, the 

• 2 

metric signature is (— , +, +, +) and R = 6(^ + ^) is the Ricci scalar, a is the scale factor and a 
dot denotes differentiation with respect to the cosmological time. 

After dropping the full derivatives with respect to time we obtain the dynamical equation for 
scalar field from variation 6S/6(p = 

(f> + 3H^ + ^R4) + eU' ((/)) = 0. (5) 

as well as the energy conservation condition from variation 6S/6g = 



£ = eh^ + £3^H^(f>'' + £3^i7(02)-+ U{<f>) - ^H^ (6) 



If we postulate the existence of a barotropic fluid in the universe the conservation condition 
reads 

-^H'^ = P<t> + Pm, (7) 

2 

H = -Y{{P4'+P^)+Pm{l+Wm)) (8) 

from which we can simply receive the energy density and the pressure of the scalar field, namely 
p^ = e^</)2 + U{(t)) + £6^i?# + £3^i^V^ (9a) 
= £^(1 - 40<^2 _ ^ £2^H4>^ - £2^(1 - e^Hcl)'^ - £3^(1 - SOi^V^ + 2^0;/' (</)). (9b) 
and in the flat FRW model the equation of state (EoS) for the scalar field is given by 

= ^. (10) 
P4> 

In the case of non-minimal coupling ^ ^ the additional variable z should be chosen in the 
parameterization of dynamics 

(U, 
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The choice of phase variables is suggested as in the case of minimal coupling by the energy conser- 
vation condition written as 

2 2 

or in terms of dynamical variables x, y, z 

^<t> = y^ + - 60^^ + 6C(x + zf) = 1-0^ (13) 

where the sense of variables x and y is preserved like in ^ = case. This condition defines in 
the phase space a domain 0^ > admissible for motion. If Vtm = then condition (jl3p define a 
2-dimensional surface in the space {x,y,z). 

Similarly the acceleration equation can be rewritten to the form 

3 

H = - — {p^s + Pea) = -^H\l + w,s), (14) 

where 

Pes Pel, + WmPm I^^P^ , o /1 c\ 

WeS = = = ^TTTT + WnMm (15) 

and after the substitution of ()14|) to (I9b[) expressed in terms of variables (x, y, z) we obtain 

2 

HI = e{l - 4e)x2 - + 2eAz) + e4^xz + eUfz^ + ^6^(1 - Qi)z^WeSf (16) 

where A = ~'^jf^- Finally after substitution of above formula into (jlSj) we can explicitly 
calculate Wcs, namely 

e(l - 4e)x2 - y2(i + 2i\z) + e4Cxz + ell^'^z'^ + w^n^ 
"^^'^ 1 - .6^(1 - 60^^ ' ^ 

where is given by equation ([13]). 

Let us start to find the dynamical system describing the evolution of our model — the non- 
minimally coupled scalar field in the Robertson- Walker background. For derivation of basic equa- 
tions we take logx (or logy) and then differentiate with respect to the cosmic time variable both 
sides of the expression. Then we obtain 

-=i-i^-,H-iRt-.m-^ (18) 

X d) H (j) (j) H 

where we have used the dynamical equation for the motion of the scalar field (j) + 'iH(j) + (,R4> + 
eU'{(p) = and R = 6(| + ^) = 6{H + 2H'^) is the Ricci scalar. The above equation reduces to 

Hx X ^y/6H^x m ^ > 
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after the reparameterization x = H-^^ we obtain 

dx „ 1, 2 H 



dlna ^ ' 2 ^ H"^ 

Finally, we obtain 



3x-12ez + e-Ay^-— 2(x + 6Cz) (20) 



dlna 2 

1 - e6e(l - Si)z^ + e{l - 4e)2;2 - y'^{l + 2^Xz) + eA^xz + wim^^l}) 



3 X + 6f z 



2 1 - e6C(l - 6C)z2 
where is determined from the constraint condition (jl3[) . 



The same method can be adopted to the variable y = ^^^j namely 



y 2 " 

and 

dy ^ _1 
dlna 2 

3 y 



y_iuw, H ^^^^ 



Xxy 



2 1 - e6C(l - 6C)z2 

And after the elimination of 0,^ equations (j2ip and (j23p can be presented in the form 



'1 -e6C(l - 6C)z2) = -3x-12^z + eiAy2(i_e6C(l-60) +e6C(l-6C)xz2 + 

+^{x + 6^z) - 6^(1 - u'^)x^ + e2C(l - 3w;™)(x + z)^ + 



dlno 



+ {l+Wrn){l-y^)) (24) 



^1-66^(1-602^) = --Ay(x(l -£6^(1- 6^)2^) +6Cy2^) -£12^(1-60^2^ + 



dlna^ .V y 2 

+^y(e(l - 60(1 - W'm)x2 + e20l - 3w;™)(x + zf + 

+{l + wm){l-y^)) (25) 

where Wm = 0, 1/3 for dust matter and radiation respectively. 

The dynamical equations (p^ and ([25|) should be completed by two additional equations to 
make the dynamical system closed, namely 

dz _ 
dlna 

and the last equation can be established from the definition of A variable 



(26) 



dA _ ,2 
dlna 



A2(r-l)x (27) 
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where T = 



U"U 
-Jjrr- 



Making following time reparameterization 



(1 - 66^(1 -60^^)^ 



d7 



(28) 



we can write complete dynamical system in the form 



X 



y 



z 



-Sec - 12^2 + £- Ay2 (1 - £6^(1 - 6^)) + £6^(1 - 6C)xz^ + 

+ ^{X + Q^Z) {e{l - 60(1 - Wm)x' + £2^(1 - ^Wm){x + zf + (1 + Wm){l 

-\\y{x{l - £6^(1 - eO^') + e^y'z) - £12^(1 - G^y^^ + 

+ \y{e{l - 60(1 - Wm)x^ + £2^(1 - ^Wm){x + zf + (1 + Wm){l - 2/')) , 

x(l-£6^(l-60-2'), 
-\\T-l)x{l-eQi{l-Qi)z^) 



y')),(29a) 



(29b) 



(29d) 



(29c) 



It is also another way to eliminate one of the variables, namely z. If we assume that V = T{X) 
then the variable z can be expressed by A according to the formula 



If we assume that r(A) = 1 — a/X^ for an arbitrary constant a, then z{X) can be integrated in the 
exact form 

.(A) = A. 
a 

Let us now make some important remarks about the general properties of the system describing 
the evolution of the non-minimally coupled scalar field on the background of the Robertson- Walker 



1. In the general case of non- vanishing barotropic matter the system is reduced to the form 
of a 3-dimensional autonomous dynamical system which can be studied by the dynamical 
systems methods. The motion of the system is restricted to the domain ^l^{x,y,z) > 

2. If effects of the barotropic matter are not considered then motion of the system in the phase 
space is restricted to the 2-dimensional surface = 1 

3. If the form of the potential function is assumed at the very beginning then the scalar field 
cosmological model is represented by a 2-dimensional dynamical system. 




(30) 



symmetry: 
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III. DYNAMICAL SYSTEM WITHOUT THE MATTER n,n = 

Let us consider the model without matter. The motion of the system takes place on the surface 
determined by equation (11). Putting $7^ = this surface is given by 

n^ = l =^ y'^ = l-e{{l-6C)x^ + 6Cix + zf). (31) 

By using this equation we can eliminate the variable y from equation ()17p and now the effective 
equation of state coefficient assumes the following form 

"^^^ ^ 1 - £6^(1 - gg)^2 {g(l-^g)^'+g^^^^ + gl2g'^'-(l + ^^^^)(l-g((l-60^' + 6g(3; + ^)'))} 

(32) 

If we assume that F = r(A) then from equation (27) we can eliminate the variable A and reduce 
the dynamical system to the two-dimensional one 

^ = _x(l-e6C(l-60z2^ +e(l-6e)(x + 6e2)x2 + 

+ (l - e((l - 60x2 + 6C(x + z)2)) {e^X{z){l - e6^z{x + z)) - 2(x + 6^^)) , (33a) 

^ = x(l- £6^(1-60^2). (33b) 
dr 

Dynamics of this system for a quadratic potential function X{z) = —2^ and the phantom scalar 



2l| 



field e = —1 have been studied in our previous work 

Another possibility is the elimination of the variable z and then the dynamical system is in the 
form 

^ = -x(l-e6ai-6e)z2(A)) +e(l-6e)(x + 6^z(A))x2 + 

+ (l - e((l - 60x^ + 6C(x + z(A))2)) (eiA(l - eQCz{X){x + z(A))) - 2(x + 6Cz{X)J^^a) 

^ = -A2(r(A) - l)x(l -£6^1 -60^^). (34b) 

For the model with $7^ = and arbitrary coupling constant we have following critical points 
at finite domain of the phase space 

• (xo = 0, Ao = : z{Xo) = 0) 

the effective equation of state parameter calculated at this point is 

WeS = -1 (35) 



the eigenvalues of the linearization matrix are 

= _l(3 + ^9 + e2a-48e), h = -^(3 - ^9 + e2a - 48^) (36) 
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-1 -0.5 □ 0.5 1 



FIG. 1: Phase portrait on the plane (x, A) in the finite domain for model without the matter f2„j — and 
r(A) = 1 — as an example. The critical point in the origin (0, 0) is of a stable focus type (see Table 

where a = /(Aq) = -A2(r(Ao) - 1) 

• (xo = 0, Ao = const. : z^(Ao) = ^) 

the effective equation of state parameter is 

^^off = ^ (37) 

the eigenvalues of the linearization matrix are 

^1 = -6?, h = 24e (38) 

The stability of the critical points depends on the sign of the real part of the eigenvalues. The 
first critical point will be always stable if only 9 + 2ae — 48^ < 0. For the second critical point 
which represents the radiation dominated universe we obtain a saddle critical point. 

IV. MODEL WITH THE MATTER 0„ 7^ 

Let us consider the general case of the model with barotropic matter and the scalar field with 
non-minimal coupling. The evolution of the model is represented by a 4-dimensional dynamical 
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TABLE I: The simplest finite critical points for model with no matter content = and arbitrary coupling 
constant ^. 





Critical point 


Wcff 


eigenvalues 


existence 


1) 


xo = 0, Ao = 


-1 


l-^ = _i(3 + V9 + s2a - 480 
^2 = -i(3- V9 + e2a-480 


z(Ao) = 


2) 


a:;o = 0, Ao = const 


1 
3 


h = -6C 
h = 24^ 


eC > 



system with polynomial right-hand sides 

nm = l-y''- 4(1 - 60a^' + 6^(x + z)^) 

and effective equation of state parameter in terms of dynamical variables reads 

1 



WeS 



1 _ £6^(1 - 60^{ - 1 + - - "'"^)^' + ^2^(1 - 3wm){x + zf+ 
+ (1 + ii;^)(l - y^) - £2^(1 - 60^' - 2^\y''z], 



Dynamical system is in the form 



dx 
d7 



dy 
dr 



dz 



-3x - 12^z + s^Xy^{l - £6^(1 - 60) + £6^(1 - QO^z^ + 



(39) 



(40) 



+-(x + 6^z) (£(1 - 60(1 - Wm)x'' + £2^(1 - 3wm)ix + zf + (1 + Wm){l - y')) ,(41a) 
-iA2/(x(l - £6^(1 - 60^') + Qiy^z) - el2i{\ - d^yz^ + 

+ ^y(£(l - 60(1 - Wm)x'' + £2^(1 - 3Wm){x + zf + (1 + Wm){l - y')) , 



— = a;(l-£6e(l-60^'), 
^ = -AV-l)x(l-£6^(l-60.^) 



(41b) 
(41c) 
(41d) 



In the special case of a minimally coupled scalar field the dynamical system is reduced to the 
following 

dx 



din a 

din a 

Az 

din a 
dA 

din a 



-3x + £^A/ + ^a;(£(l - Wm)x^ + (1 + Wm){l - y^)) , (42a) 

-^Xxy + ^y(£(l - Wm)x^ + (1 + w;^)(l - y^)), (42b) 

X, (42c) 

-A2(r-l)x. (42d) 
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In what follows we will assume that T = T(\) and that for any critical point exists Aq such that 

/(Ao) = A2(r(Ao) - 1) = -a (43) 

and 

^|ao = /'(Ao) = const (44) 

are finite. 

For non-minimally coupled scalar field in the energy phase space we have following critical 
points: 

• xq = 0, uq = 0, and Aq = const, where Aq is an arbitrary constant such that 2;(Ao) = 0: 
the effective equation of state parameter calculated at this point is 

WeS = Wm (45) 

the eigenvalues of the linearization matrix are 



4 V y 3 [l-WmY 

h = ^{l+wm), (46a) 



4^ '"^V Y 3^(l-u;^)2 

For barotropic fluid equation of state parameter Wm > — 1 the eigenvalue I2 is always positive 
and this critical point corresponds to an unstable focus when the eigenvalues li and ^3 are 
complex numbers or to a saddle when li and ^3 are purely real and negative. In opposite 
case, when Wm < — 1 this point correspond to a sink (i.e. a focus or a stable node, depending 
on the value of square root in li and Z3) because the real parts of the eigenvalues are always 
negative. 

{xo = 0, y2 = 1, Ao = : z(Ao) = 0) 

the effective equation of state parameter at this point is 

Weff = -1 (47) 

the simplest form of the function r(A) which fulfills both conditions and (HH) is 

r(A) = 1 - (48) 
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TABLE II: The simplest finite critical points for dynamical system describing model with matter flm 7^ 
and arbitrary coupling constant ^. 





Critical point 


WeS 


eigenvalues 


existence 


1) 


xo — 0,yo — 0, Ao = const 




h = -1(1 - w^)a + ^1- f ^(Hjf^) 

h = 1(1 + w„0 

V 


^(Ao) = 


2) 


xo = 0, = 1, Ao = 


-1 


li = -i(3 + V9 + £2a-480 

h = -3(1 + w„i) 
/3 = -1(3- V9 + e2a -480 


z(Ao) = 


3) 


xq — 0,yo ~ 0, Ao = const 


1 

3 


^1 = -ee 

^2 = 12^ 
^3 = 6^(1 - 3Wm) 


> 



the eigenvalues of the linearization matrix are 

h = -^(3 + V9 + e2a-480, ^2 = -3(1 + w;™), ^3 = -^(3 - ^9 + e2a - 48^) (49) 
the character of this critical point depends on the value of a and ^. 

• (xo = 0, yo = 0, Ao = const : z'^{Xq) = ^) the effective equation of state parameter is 

1 . , 

WeS = 2 (50) 

and the eigenvalues of linearization matrix are 

h = -6^, 12 = 12^, l3 = 6a'^-3wm) (51) 

At this critical point the effective gravitational constant changes sign an any FRW model 
becomes unstable with respect to arbitrary small anisotropic or inhomogeneous perturbations 
and a curvature singularity forms in this point as it was shown first by Starobinsky [48.] . 

V. CONCLUSION 

In this paper we extended the analysis of the dynamics of the FRW model with a minimal cou- 
pling to gravity scalar field (both canonical and phantom) to the case of a non-minimal coupling. 
We showed that in the general case the dynamics can be represented by a four-dimensional dynam- 
ical system. However, with the assumption that the form of the potential parameter is T = r(A), 
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FIG. 2: Three dimensional phase portrait of the system (|^T|l . For illustration we put r(A) = 1 — ^ and 
the same values of ^ and a as in Fig. [TJ The figure illustrates a quintessence multiple scenario with a stable 
focus type critical point sf as a final state which is the deSitter attractor {wcH — —1) and a saddle type 
critical point s (wcs =1/3) and an unstable focus uf (t«cff = Wm) as intermediate states. 

the corresponding system can be reduced to the form of a 3-dimensional dynamical system. Using 
the dynamical system methods we analyzed critical points appearing at a finite domain of the 
phase space and trajectories in their neighborhood. These trajectories can be obtained from the 
linearization of the system around a critical point of the type being determined by eigenvalues of 
the linearization matrix. 

We found some additional points in comparison to the case of the minimal coupling. There are 
in principle at most four families of critical points which character (stability) depend on the value 
of the non-minimal coupling constant. In the special case of ^ = this system was investigated 
by Zhou [46] but we found that the critical points which were established by them lie rather at 
infinity {(f) = oo). These critical points are out of interest because (p = oo may lead to singularities. 
However, we use the same methodology as Zhou and instead of postulating quintessence potential 
directly the relation between T and A is proposed. 

For our dynamical analysis it would be useful to distinguish two cases: with matter 
and without matter 0^ = 0. The former is corresponding to a 3-dimensional dynamical system; 
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FIG. 3: Plot of the evolution of Wcs (the relation (|40)) ) for the non-minimally coupled canonical scalar field 
and positive coupling constant ^. The probing trajectory used to plot this relation starts its evolution at 
To = near a saddle type critical point {wcs = 1/3) and then approaches an unstable focus critical point 
Weft — = — 1/3 and next escapes to the stable de Sitter critical point with Wcn = — 1. Note the existence 
of a short time interval during which Wcff — ^ ■ 

in the latter one variable is eliminated due to the constraint condition Q{x, y, z) = 0. For the case 
with matter we have three families of the critical points 

1. (0, 0, 0) — the barotropic fluid dominated universe which is unstable if Wm > —1 and stable 
if Wm < -1- 

2. (0, ±1,0) — the de Sitter attractor (repellor) which is unstable if Wm < —1 and stable if 

Wm > -1. 

3. (0, 0, ±-7=) — it is a radiation dominated universe of the saddle type; here are the conditions 
C / and > 0. 

We can see that the barotropic fluid dominated universe is stable (unstable) when the de Sitter is 
unstable (stable). 

From the dynamical analysis of the 3-dimensional dynamical system (see Fig. [2]) it is found a 
generic quintessence evolution scenario (in the sense that it is realized by a very wide range of 
initial conditions). In this scenario the final state is represented by the de Sitter attractor with 
Wed = — 1 (see Fig. [3|) or the Einstein-de Sitter universe with Wm > — 1 and Wm < —I respectively. 
Therefore all roads lead to the quintessence model and give the current acceleration. 
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Let us consider the trajectory starting from the point {xo,yo,zo) then it goes to a close neigh- 
borhood of a saddle point (a radiation dominated universe) then approaches to stable critical point, 
i.e. for the Wm < — 1 matter dominated universe or Wm > — 1 de Sitter state. In the latter the 
trajectory goes close to a barotropic fluid dominated universe (unstable point) before launch to 
the de Sitter state. 

For the case without the matter 0,^ = we have found the new generic quintessence scenario 
which appears only for non-vanishing coupling constant. In this scenario in the phase variables 
(see Fig. [1]) trajectories spend long time in the neighborhood of the saddle and then escape to the 
de Sitter attractor which leads to accelerated expansion. 
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